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We apply the renormalisation group (RG) to analyse scattering by short-range 
forces in systems with coupled channels. For two <S-wave channels, we find three 
CN \ fixed points, corresponding to systems with zero, one or two bound or virtual states 

f-H | at threshold. We use the RG to determine the power countings for the resulting 

effective field theories. In the case of a single low-energy state, the resulting theory 

O ; 

takes the form of an effective-range expansion in the strongly interacting channel. 

Ph . 

We also extend the analysis to include the effects of the Coulomb interaction between 

(N 

£> ' charged particles. The approach is then applied to the coupled p + Li and n + 'Be 

t> : 

0^ . channels which couple to a J p = 2 state of 8 Be very close to the n + 7 Be threshold. 

o 

At next-to-leading order, we are able to get a good description of the p + 7 Li phase 
shift and the Be(n,p) Li cross section using four parameters. Fits at one order 
higher are similarly good but the available data are not sufficient to determine all 
five parameters uniquely. 

I. INTRODUCTION 

Effective field theories (EFTs) now provide a widely used tool in studies of a variety of 
systems in nuclear, particle, and atomic physics (for reviews, see Refs. 1 3)] ) . The idea 
behind them is that, in order to describe processes at sufficiently small energies, one does 
not need to treat explicitly the underlying short- distance physics. Provided there is good 
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separation between the energy scales of interest and those of the short- distance physics, 
one can work with the appropriate low-energy degrees of freedom and expand the resulting 
theory in powers of ratios of low-energy to high-energy scales. A key ingredient of any such 
theory is thus the "power counting" used to organise this expansion of its interactions. 

Some of the main applications of these ideas have been to few-nucleon systems on momen- 
tum scales of the order of the pion mass. These make use of Chiral Perturbation Theory to 
determine the long-range pion-exchange forces between nucleons, and replace the unresolved 
short-range forces by contact interactions. At lower energies, even pion-exchange forces are 
not resolved and one can work instead with a simpler, pionless EFT, involving only contact 
interactions [4]. In nucleon-nucleon S-waves, the deuteron bound state and the virtual 
state have the effect of enhancing the low-energy wave functions at short distances. The 
resulting power-counting in these channels shows a strong promotion of the leading contact 
interactions, which must be treated nonperturbatively. 

More recently, a similar approach, also based on contact interactions, has been applied 
to larger nuclei which have low-energy states or resonances with a clustered structure. This 
theory, known as "halo EFT" , has been applied to various weakly bound nuclear systems 5^ 
9[. A number of these applications are to systems of astrophysical importance, and some of 
the most interesting of these involve reactions whose cross sections are enhanced by resonant 
or virtual states lying very close to a threshold. EFTs that incorporate such states have 
been developed in Refs. lOl 111 ]. 

n 

Theories describing coupled channels, such as those developed by Cohen et al. [1Q|] are 
needed for applications to reactions. Here we use the renormalisation group (RG) to de- 
termine the possible power countings for systems with two scattering channels. This is an 



extension of the Wilsonian approach 



12] , introduced in Ref . 13 1 to study single- channel sys- 



tems with short-range interactions. We first determine the fixed points of the RG, and then 
use the linearised RG equation to analyse the scale dependences of perturbations around 
these points and hence find the power counting. 

In the case of two coupled channels, we find three fixed points. One of these is just the 
trivial one, describing a non-interacting system. This is the appropriate starting point for an 
EFT describing a weakly interacting system, with no enhancements of the scattering near 
thresholds. A second has strong interactions in both channels, corresponding to two virtual 
states lying at the higher threshold. This corresponds to the case considered in Ref. 101 ]. 
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where the matrices of energy-independent interactions are promoted compared to simple 
dimensional power counting. 

Lastly there is a fixed point with a single virtual state at the threshold. This is the 
one of most physical relevance since it requires "fine-tuning" of only one quantity. This 
is in contrast to EFT studied in Ref. 10] where two parameters, the eigenvalues of the 
leading matrix interaction, must be fine-tuned. Both that EFT and the one studied here 



correspond to coupled-channel extensions of the effective-range expansion 

H4H16|. They 

provide potential alternatives to the i?-matrix approach {l^] that is widely used to analyse 
reactions with resonant or virtual states fl8 |. 

As well as being applicable to low-energy resonances in nuclear physics, the same coupled- 
channel EFT can describe near-threshold states in other contexts, such as quarkonium states 

n 

in hadronic physics. One particularly interesting example is the X(3872) [19], which lies 
very close to the D°D*° threshold. This was recently studied by Hanhart et al. using a 
phenomenological coupled-channel approach 20]. It would be very interesting to examine 
it within the EFT framework developed here. 

In general, at least one of the channels in a nuclear reaction involves two charged particles. 
As well as the short-ranged nuclear forces, we therefore need to treat the Coulomb inter- 
action. EFT and RG techniques have been developed to analyse scattering in the presence 



of a Coulomb potential 



2lN25|. These show that the 1/r singularity is not strong enough 



to alter basic power counting, which is still that of an effective-range expansion. However 
there is a new low-energy scale k, the inverse Bohr radius. The interaction also leads to 
a logarithmic divergence that needs to be renormalised by a counterterm linear in k. This 
logarithmic behaviour makes it impossible to disentangle the purely strong-interaction scat- 
tering length from the scatttering data. In the present work, we extend these treatments to 
coupled-channel systems, concentrating on the case where only one of the channels consists 
of two charged particles. 

As a practical illustration of the resulting EFT, we apply it to the system of coupled 
p + 7 Li and n + 7 Be channels. The reaction 7 Be(n,p) 7 Li is an important one in the con- 
text of primordial nucleosynthesis of 7 Li. However, it is one that has been well studied 

3o| and so it is not expected to provide an explana- 



experimentally over many years 26 



tion for the observed abundance of 7 Li 



3jj. The astrophysical importance of this reaction 



is a consequence of its enhancement by a virtual 2 state of 8 Be lying very close to the 
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n + Be threshold. This makes the system an ideal one to test our EFT approach, as it has 
a single low-energy state which is coupled strongly to both physical channels. We calculate 
scattering at next-to-next-to-leading order (NNLO) in our power couning and determine 
the corresponding low-energy parameters from the available experimental data. We also 
examine the nature of 8 Be pole within our approach. 

The structure of our paper is as follows. In Sec. [Til we use the RG to analyse scattering 
in systems with two coupled channels. The extension of this to systems with Coulomb 
interactions is outlined in the Appendix. In Sec. IIHI we apply the resulting EFT to the 
example of coupled p + 7 Li and n + 7 Be channels. We close with some conclusions in Sec. IIVI 



II. RG ANALYSIS 

A. Short-range forces 

Here we use the RG to analyse scattering in a two-channel system with only short-range 
forces. The Lippmann-Schwinger equation for the T-matrix can be written schematically as 

T = V + VGT, (1) 

where 

(2) 




in terms of the free Green's functions for the individual channels. These can be written in 
momentum space as 



Gi = 2M, / ^ , (3) 



d 3 q 

2ir 3 pj — q 2 + i e ' 

where Mj is the reduced mass in channel % and pi is the on-shell momentum. Taking zero 
energy to be the lower threshold and A to be the threshold energy of the second channel, 
we have 

pi = y/2M 1 E = p, (4) 



p 2 = ^2M 2 (E-A) = ^j^-(p 2 -5 2 ), (5) 
where the momentum scale 5 associated with the difference between the thresholds is 

5 = a/2MiA. (6) 
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Specialising to the case of two S'-wave channels, we can take the short-range potential to 
consist simply of 5-functions with energy- dependent coefficients. This is because terms that 
depend on the off-shell momenta (that is, derivatives of 5-functions) are of the same or higher 
order as the corresponding energy-dependent ones and so are not essential for a description 



of the on-shell scattering 13|, |24], [25] . This leaves the system with two low-energy scales: p 
and 5. These provide the expansion parameters of our EFT. 

For S'-wave channels with contact interactions, the coupled integral equations for T reduce 
to algebraic equations for the coefficients of the (^-functions, 

T(p, 5) = V(p, 5) + V(p, 5)J(p, 6)T {P , 6), (7) 

where the diagonal matrix J(p, 5) consists of the loop integrals 

J t (p,6) = 2M t [ ff 2 1 ■ (8) 
J 2ti a pf — q z + i e 



These loop integrals are divergent and so need to be regularised in some way. In Ref. [131 ]. 
a sharp momentum cutoff was used. The resulting expressions for the potentials are a 
little cumbersome and so here we use dimensional regularisation with the power divergence 
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33| . As discussed in Ref. |l3| , 



subtraction scheme intoduced by Kaplan, Savage and Wise 
the results for "universal" quantities — RG fixed points and power countings — are the same 
as those obtained with a cutoff. 

In the absence of Coulomb interactions, there are only linear divergences and subtracting 
these at the scale p gives 

Ji{p,S,n) = -—^{n + ipi). (9) 

The subtraction scale /i is arbitary and so the physical scattering amplitudes in T should be 
independent of it. To cancel the \i dependence of the loop integrals in Eq. (jSJ), the potential 
matrix V(p, S, p) must run with p according to 

To determine the possible power countings, we need to put Eq. ( FlO|) into the form of a 
standard RG equation. We do this by defining dimensionless variables corresponding to the 
low-energy scales, 

P = p/p, 5 = 5/p } (11) 
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and a rescaled potential, 
where 



V^M^VM 1 ' 2 , (12) 

27T 



■ -Mi , 

M = . (13) 

M 2 

The evolution equation then becomes 

<9/i <9p 55 

The dimensionless form of this equation allows us to look for fixed-point solutions which 
describe scale-free physical systems. We then expand the potential around one of these points 
and find perturbations that scale with definite powers of the regulator scale \x. Because of 
the rescaling, this power of [i counts the net power of low-energy scales in the corresponding 



term in the potential and so gives the order of that term in the power counting [13j . More 
precisely, a term that scales as fi u is of order Q v ~ 1 1 where Q denotes a generic low-energy 
scale. 

One obvious fixed point is the trivial solution to Eq. (TH|) : V = 0. The EFT based on 
the expansion around this point decribes systems that interact weakly at low energies. The 
leading (energy-independent) term in the potential scales as /i 1 and so is of order Q°. Each 
power of the energy (p 2 ) and splitting (5 2 ) increases the order by Q 2 as expected from naive 



dimensional analysis (or "Weinberg power counting") 34j. |35|. 

A second, nontrivial fixed point can easily be constructed by analogy with that in the 
single- channel case. If we express the RG equation in terms of V -1 , then it takes the simpler, 
linear form, 

A* ^V- 1 = p ^V- 1 + 8 ^V- 1 - V- 1 - I, (15) 
dfi ^ dp 05 

where I is the 2x2 identity matrix. This has the //-independent solution 

V = -I, (16) 

which is just the 2 x 2-matrix equivalent of the fixed-point that corresponds to the effective- 
range expansion for a single channel jl3]. It describes a system with bound states at threshold 
in both channels. 
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Since Eq. ([15]) is linear, exact solutions can be constructed and expressed in terms of 
perturbations that scale with definite powers of fi: 

V(p, 6, fi)' 1 = -I - C nm ^ n+2m ~ 1 p 2n 5 2m . (17) 

n,m 

This shows that the fixed point has two unstable directions, given by the two eigenvectors 
of Coo- The expansion is the appropriate one for systems where both eigenvalues of Coo, 
denoted by Ci^, are unnaturally small, corresponding to two low-energy bound or virtual 
states. This is the counting considered by Cohen et al. 1(3], who treated all elements of the 
energy-independent matrix Coo as leading-order in their expansion. In this case we can keep 
/i > c\ t 2 so that the constant —1 provides the leading term in V -1 . Taking the inverse, we 
get the potential in the form 

V(p, 6, fj) = -I + Y, C nm /i 2 ^ 2 ™" 1 p 2n S 2m + ---, (18) 

n,m 

which shows that all terms are enhanced by two orders compared to Weinberg power count- 
ing. 

Of more interest is a second nontrivial fixed point, which has the single-term separable 
structure 

V = uK,(p,5)ut, (19) 



where 




u= Y \. (20) 



Its strength V (p,S) satisfies 



and hence is just 



V (P, <*) = -!■ (22) 

This potential generates a single zero-energy bound state whose overlaps with the two scat- 
tering channels are given by the components of u. In the orthogonal combination of channels, 
there is no interaction. 

Physical systems with a single low-energy bound or virtual state can be described by 
potentials close to this fixed point. Expanding V around it and linearising the RG equation, 
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one can show that interactions in the channel u are enhanced by two orders, as in Eq. ([IB"]) 
above. In contrast, those in the orthogonal channel, specified by 



(23) 



have natural coefficients and can be organised according to Weinberg's power counting. 
Terms in the potential that couple the u and v channels are enhanced by one order. 

In practice, an easier way to construct the expansion of a potential around this point is to 
start from the general solution to the RG equation, Eq. ( ITT)) . If only one of the eigenvalues, 
Ci, of Coo is unnaturally small, we should choose our renormalisation scale fi such that 
c\ < fi <C C2- In this case, the term C2//X becomes the dominant one in the expansion of 
V -1 . As a result, the expansion of V has a different structure from the case just discussed, 
Eq. ffTB"]) . In particular, the expansions in the channels corresponding to the two eigenvectors 
of Coo have different power countings ^39 ] . 

To find the corresponding expansion of the whole potential, we need to invert Eq. ( fl7j) . 
It is convenient to introduce four matrices P a , defined by 

P u = uu\ P v = vv t , P m = uv t , P] n = vu t , (24) 

where u and v denote the eigenvectors of Coo corresponding to the eigenvalues c\ and C2 
respectively. In terms of these, we can write the general solution for the inverse potential, 
Eq. (|T7|) . in the form 



V(p, 8, /i)- 1 = f n (p, 8, /i)-^ + / v (p, 8, /U )- 1 P V + f m (p, 8, fiy'Prn + /* (p, 5, ri-'Pl, (25) 
where the functions f a (p, 8, /x) _1 are defined by 

n,m 

UPXhY 1 = "I - J2°nl /i 2 "^- 1 ^"^, (27) 

n,m 

UP,5^)- 1 = -^c^^+ 2 - 1 p 2 "5 2 - (28) 

and Cqq = Ci, Cqq = c 2 . Note that = 0. Provided no other channels are open, this 
potential should be Hermitian, with real. 
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Inverting Eq. (|25|) we get the potential, 



V \_f v Pu + f n Pv /m Pm / m ^m] [/u /v fia fm ] 



(29) 



For the situation of interest, we need to expand this assuming that c\ <C fi <C c 2 . In the 
channel corresponding to the small eigenvalue, this gives the diagonal interaction 

= [/u ~~ fv fm fm ] 



1 _\^ r ( u ) ,.2n+2m~l ^2n pm j_ „(m)^2 , „(m) £ 2 2 



C 2 



+ 



(30) 



where terms suppressed by additional powers of /i/c 2 or c\j\i have been omitted. The 
dominant terms multiplying each product of powers of p and 5 have the same power counting 
as in the effective-range expansion, that is, enhanced by two orders over simple dimensional 
analysis. 

In contrast, the diagonal interaction in the other channel has the expansion 



-l r*— in — i 

' m J 



r ( m V , (m) ?2 
c 10 V < c 01 



+ 



— [f v f u f m /] 

Yl C ni ^ n+2m - 1 5 2m ~ 1 + V 2 

n,m 

£ cffi /^+ 2 ™p 2 « 5 2 ™ + A* (-1 + /x 2 jc^p 2 + c^ 2 

2 nm " 



+ 



n,m 

ii _ 
C 2 C2 



(31) 



,2n+2m+l ~2n r2m 



where the prime indicates the omission of the term with n = m = 0. This is just the standard 
Weinberg power counting, as expected for a channel with no low-energy bound or virtual 
state. The renormalisation in this channel is done perturbatively; the terms (suppressed in 
the equation above) necessary for this purpose do not affect the power counting. Finally, 
the off-diagonal interactions can be expanded as 

V =-f~ 1 f\f~ 1 -f /■ — 1 -1 

'm J m J v U u J v J ra J m J 
_J_ \ ^ 2n+2m ~2n i2m , 

(32) 

/ , C nm h 1 P T ' " " ■ 

C 2 



n,m 



The terms in these are enhanced by one power compared to simple dimensional analysis. 



10 



In the corresponding unsealed potential, the RG eigenvalues translate to powers of low- 
energy scales, generically denoted by Q, a term that scales as \i v being of order Q d with 
d — v — 1 [l3|. The small scales here are p and 5, and also the subtraction scale fi, which can 
be regarded as the highest acceptable momentum scale in our EFT. However, terms with 
different powers of 5 cannot be distinguished in practice and so they can be grouped together 
to form observables, such as the scattering length in the strongly-interacting channel, 



1 

a i 



-ci + 0{5 2 



(33) 



The unsealed inverse potential corresponding to Eq. ( Tl7l) can then be expressed in the 
form 



1 



V- 1 = — — M 1 / 2 
2vr 



I + ^)p„ + L_ljp v + ^ ( P m + P y + ... 



M 1 / 2 , 
(34) 

where we have expanded in powers of energy around the higher threshold, with p' = p 2 , and 
kept terms to NNLO (order Q 3 ) in the power counting defined above. In this counting the 
large scattering length ai is of order while all other coefficients are of natural size, that 
is, of order Q°. We have assumed that no other channels are open and hence we have taken 
the mixing parameter r\ to be real. 

The Lippmann-Schwinger equation is conveniently written in the form 



T- 1 = V 



-i 



J. 



(35) 



where the terms proportional to [i in the loop integrals can immediately be seen to cancel 
against the corresponding terms in V -1 . The resulting expression for T _1 can then be 
inverted and expanded in the same manner as used above to construct the potential V. The 
resulting expression for the scattering matrix is, to NNLO, 

-l 



T NNLO = _ 2vr M- 1 / 2 . 



1 



r 



p' 2 - iPu - a 2 ([1 - ia 2 p v ] ^ - irip /2 p m ) 



«2 



+ a 2 



1 - ia 2 p v + a 2 p m ( -— - ip u 



a i 



- Wu 



[1 + ia 2 p v ] + r -^p' 2 - a 2 pl 



a i 



-1 



Pi) >M" 1/2 , 



(36) 



11 



where we defined the momentum variables 

2/2 2/2 / f \ / \ 

p u = P cos + p sin 0, p v = p sin <fi + p cos 0, p m = (p — p ) sin cos 0. (37) 

Here we have chosen to leave the diagonal amplitude in the strongly-interacting channel, 
as well as parts of the amplitude in the mixing channel, in the form of an effective-range 
expansion but, if desired, these terms could also be expanded to order Q 3 . 



B. Short-range plus Coulomb 

With some modifications, the analysis above can also be applied in the presence of the 
Coulomb interaction in either or both of the channels. This introduces two new low-energy 
scales, the inverse Bohr radii for the two channels, K± t 2- The Lippmann-Schwinger equation 
takes the form 

T c (p, 5, «i i2 ) = V(p, 5, «i t2 , n) + V(p, 5, «i j2 , fi)J(p, 5, n ly2 , fi)T c (p, 5, «i >2 ). (38) 
The individual loop integrals in J(p, 5, Ki$, p) are given by 



^ S (P, <5, «,/*) = - ^ (V - 2k, (in ^ + 1 - ? 7 4 + 2 Ki 



(39) 



where rji = Ki/pi, h(z) = Re^iz) — In z (ip denoting the logarithmic derivative of 
Euler's gamma function), C|. = 27r^j/(exp 2irr]i — 1) are the Sommerfeld factors, and 
7£ = —0.5772 ... is the Euler-Mascheroni constant. To get this, we have used PDS for the 
linear divergence and minimal subtraction for the logarithmic one induced by the Coulomb 
potential. 

These logarithmic divergences in the basic loop integrals generate logarithmic depen- 
dences of the potential on /i. The resulting general solution to the RG equation can be 
written 

VOU&^r 1 = -I + k2\n^-J2 G nml ii 2n+2m ^ l - x p 2n l 2m k\ (40) 

n,m,l 

where A is an arbitrary energy scale, and k, is the diagonal matrix of the (rescaled) inverse 
Bohr radii. Note that Cooi must depend on A so as to cancel the explicit dependence on A 
in the second term. 

As in the single-channel case, the only true fixed-point solution is the trivial one V = 



22 



24 



25] . This is because both of the nontrivial fixed points found above become unstable 
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as a result of the logarithmic evolution induced by the Coulomb potential. Nonetheless, 
we can still construct solutions that describe systems with either one or two bound states 
that are close to threshold on the energy scales of interest. Again, it is the expansion of 
the potential around the solution with a single such state that is of more interest since it 
requires only one "fine tuning", corresponding to one unnaturally small eigenvalue of Cooo- 
The power counting for this expansion is similar to that derived above. A minor modification 



is the presence of the short-distance interactions proportional to powers of k, althoug 



i these 



21 



23|. 



cannot be disentangled phenomenologically from purely strong-interaction terms 

The corresponding scattering matrix has a very similar form to that in the purely short- 
range case. Expanded to NNLO, it can be written 

~i -i 



-.NNLO 



-2nM- 1/2 - 



1 

a i 



yP' 2 - Ju(p) + a 2 (j'm(p) [1 - a 2j v (p)} + np' 2 j m (p)) 



Pi 



a 2 



+ a 2 



«2jv(p) - Otfmip) 



1 



kip) 



kip) 



H V 



1 

1 

a i 



- kip) 
kip) 



[1 + a 2 j v {p)) + ^p' 2 + a 2j 2 m ip) 



-i 




(41) 



where we have used the fact that the linear and logarithmic dependences on \x cancel in the 
Lippmann-Schwinger equation Eq. fl35|) . This allows us to express the result in terms of the 
finite parts of the loop integrals in the presence of the Coulomb potential, 

jiip) = iu(2h{r] i )+i^y (42) 

where i — 1,2 denote the physical channel, and the relevant combinations of ji(p) are defined 
analogously to the purely strong-interaction case: 

kip) =hip) cos 2 <p + j 2 (p) sin 2 0, 

k ip) =h ip) sin 2 <p + h ip) cos 2 0, (43) 

kip) = (ii(p) - hip)) sin cos 0. 
For the NLO expression, the terms with r\ should be omitted, as well as the terms propor- 
tional to a 2 in the expansion, and the term with r in the mixing channel; for the LO one, 
the terms with a 2 and the remaining term with r should also be left out. 
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III. p+ 7 Li«n + 7 Be COUPLED CHANNELS 

A. EFT description of the data 

To illustrate the application of an EFT based on the power counting outlined above, we 
examine here the p + 7 Li -h- n + 7 Be coupled-channel system. The n + 7 Be threshold lies 
at A = 1.6442 MeV and there is a J p = 2~ excited state of 8 Be within a few keV of this 



threshold 



28 



i elastic 
for the 



29] . This state manifests itself as a prominent peak in the p + 7 
scattering cross sections [stJ and a very large cross section, a ~ 38.4 x 10 3 b {28 
reaction 7 Be(n,p) 7 Li at low (thermal) energies. Together with the absence of any unusual 
features in the 5 S 2 phase shift below the n + 7 Be threshold [27], this makes the system an 
ideal candidate to be studied with our approach. 

The system has two important low-energy scales: the momentum scale corresponding to 
the difference between the thresholds, 5 = 51.95 MeV, and the inverse Bohr radius for the 
p+ 7 Li channel, k\ = k — 17.96 MeV. These, together with the on-shell relative momentum, 
provide the expansion parameters of the EFT. 

The scattering matrix for this system is given by Eq. (I4"T|) . with the simplification that 
J2G 5 ') = ip' since the neutron in the second channel is uncharged. 

We now confront our EFT with the available data in the p + 7 Li and n + 7 Be channels. 
This data consists of the 5 S 2 phase shift, Sq, in 7 Li + p elastic scattering, where results of 
partial- wave analyses and R- matrix fits exist |27J , and the cross section, u„ of the charge 



exchange 7 Be(n,pV 7 Li reaction, from high quality measurements of Ref. 



and the R- 



matrix fit of Ref. 29J]. These quantities are related to the corresponding elements of the T 
matrix via 



pe 2l5 ° ~ 1 M x 2 
2ip " ~2^ TuC ^ 



-^T llC l (44) 



4tt(2J + 1) p M,M 2 l 2 
anp ~ (2 Sl + l)(2 S2 + l)^^^ | -^ T211 ' (45) 

where J = 2 is the total momentum in the partial wave considered, Si = 3/2, s 2 = 1/2 
are the spins of initial particles, rj = n/p, and the Sommerfeld factors arise due to the wave 
functions of the corresponding Coulomb scattering states. Here, p is the measure of elasticity 
in p + 7 Li channel, equal to one below n + 7 Be threshold. 

Before we proceed to the description of our results, a few words on our fitting procedure 
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are in order. We fit NLO and NNLO results to the phase shift and the reduced (n,p) 
cross section, <7^i d = a np \/E — A, while the LO parameters are extracted from the threshold 
values of the phase shift and the reduced cross section, rather than fitted over the whole 
region (such LO fits give similar values of parameters, and reproduce the observables with 
comparable quality). The results of phase analysis of Ref. [27] do not provide an uncertainty 
for the phase shift. We therefore assign an estimated uniform uncertainty of ±5° to the data 
points. This estimate is based on the scatter of the data points and is used in the chi-square 
fitting function. We do not aim to reproduce the features of 6q at energies above the neutron 
threshold, for the reasons that will be discussed below. As to the cross section, we take the 
experimental data and uncertainties on the total cross section from Ref. j^sj], appropriately 
rescaled to obtain the reduced cross section. Except at LO, we performed fits of two kinds - 
in one case, we fitted the expanded T matrix, as in Eq. f HTj) . whereas the other fit employed 
the exact T matrix, obtained by inverting the Lippmann-Schwinger equation with V -1 at 
a given order. This form preserves unitarity exactly. In the fits of the expanded T matrix, 
we also investigated the role of constraining them to satisfy unitarity below the threshold, 
which was done by imposing a constraint that p be close to one. At each order we fit all the 
parameters that enter V -1 at that order. 

Before presenting our results, we should comment on some of the numerical details of the 
different fits. At NLO, all the versions of our fits resulted in very similar values for the low- 
energy parameters, except for 02, which was very small in all cases, although its values could 
differ by a factor of two. At NNLO, we encountered two kinds of numerical difficulties. The 
first was that some fits, namely those at NNLO with the expanded T matrix but without 
the unitarity constraint, resulted in far too high values of the reduced cross section at higher 
energies, E — A > 0.1 MeV, despite still describing the available data well. In this region 
there are no direct experimental data on the cross section for this partial wave, as a result of 
other partial waves, especially J p = 3 + , becoming important. We therefore used the results 
of the i?-matrix fit of Ref. [29( as reference points at these energies, in order to filter out 
fits that led to unnaturally large value of the cross section at higher energies. The second 
difficulty was related to the fact that already at NLO we reproduce the available data quite 
well, and hence adding a new parameter, r 1; results in a very broad minimum of the fitting 
function at NNLO. This indicates that the available data is not enough in order to constrain 
the scattering parameters beyond NLO, and, as a result, the variations of the scattering 
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parameters between the different versions of fits grow significantly compared to NLO. 

In Table |T] we show the values of the scattering parameters that resulted from our fits, 
using the expanded T matrix without the unitarity constraint. At LO, taking the value 
of 5 5*2 phase shift at the neutron threshold 83.3° and the value of the reduced (n,p) cross 
section 5.75 b MeV 1 / 2 gives 

ai = -17.76 fm, (46) 
= 46.63°. (47) 

One can see from Tabled that the low-energy parameters a 2 , r and r\ are of natural size for 
an effective theory with an underlying length scale ~ 2 fm. The values of a\ and change 
slightly when one goes to NLO and NNLO. Our results for the phase shift and the reduced 
(n,p) cross section are shown in Fig. [TJ Surprisingly, the LO description of a^p is already 
very good up to energies ~ 0.5 MeV above neutron threshold. In fact, all three curves, LO, 
NLO, and NNLO, differ only slightly and only at higher energies (even though the p + 7 Li 
phase shift 5q at these energies is not well described at LO). Furthermore, we observed that 
fitting only the threshold value of a™p without paying attention to the details of its energy 
dependence, produces, as a rule, fits of good quality. In contrast to this, the behaviour 
of the phase shift above the neutron threshold is very unstable at LO — slight changes of 
threshold values of 5 and can lead to very different behaviour of 5 . This is due to the 
closeness of 5q to 90° at threshold and a strong inelasticity in the proton channel just above 
the neutron threshold. This instability, however, was not observed in NLO and NNLO fits. 

The phase shift below the neutron threshold is reproduced quite nicely already at NLO, 
and its convergence can be seen from Fig.[U At the same time, the slope of 5q above threshold 
is not reproduced very well. We believe this to be due to the inconsistency between the data 



on the reaction 7 Be(n,p) 7 Li of Ref. 28] used in our fits and the (older) data on the inverse 



reaction, 7 Li(p, n) 7 Be, measured in Ref. [26j and used in the phase-shift analysis of Ref. [271 ]. 



The data of Ref. 26j are significantly higher than those of Ref. [28|], as also illustrated in 
Fig. [U This means that the inelasticity in the p + 7 Li channel above the neutron threshold 
was larger in that analysis than in our fits, which could explain the difference between the 
slopes. This inconsistency between the phase shifts above threshold and the data on a np led 
us to discard the phase shifts above threshold from our fits, as mentioned above. On the 
other hand, since the inelasticities inp+ 7 Li scattering are small below the neutron threshold, 
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FIG. 1: Upper panel: phase shift at different orders, fits of expanded T matrix without unitarity 
constraint. Dash-dotted curve — LO, dashed curve — NLO, solid curve — NNLO. Data: diamonds 
- phase analysis from Ref. [23] (digitized; the data uncertainties are as described in the text). 
The dashed and the solid curves are on top of each other below the neutron threshold. Note that 
only the phase shift data below the threshold are included in the fits, as explained in the text. 
Lower panel: reduced cross section of the reaction 7 Be(n,p) Li. The curves are as in the upper 



panel. Data points: circles — data from Ref. 



28j , triangles — derived from the near-threshold data 



(digitized) of Ref. 26] on the crossed reaction, 7 Li(p, n) 7 Be, hollow and filled triangles correspond 
to targets a and b of that reference, in order. 
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Order 


a\ [fm] 





a 2 [fm] 


r [fm] 


r\ [fm] 


LO 


-17.76 


46.63° 








NLO 


-19.37 


51.82° 


-1.96 


3.79 




NNLO 


-19.11 


50.45° 


-1.07 


2.58 


-5.42 



TABLE I: Scattering parameters at different orders, resulting from the fits of the expanded T 
matrix, Eq. (|4ip . without the unitarity constraint, as described in the text. 

we have no reason to question the reliability of the results of the phase-shift analysis at these 
energies. 



B. The structure of the 2 8 Be state near the neutron threshold 

Here we would like to discuss the structure of the 2~ 8 Be state near the neutron threshold. 
In our calculation, it appears as a pole in the S matrix located right above the threshold, at 
E = E r - ir/2 = 1.71 - i0.06 MeV, giving the total width T = 0.12 MeV. These numbers 
correspond to the NNLO parameters from Table [H Note that due to the expansion, the 
different channels (u, v, and the mixing channel) of the T matrix, Eq. f l41~]) . have poles at 
slightly different positions, which would not be the case for the exact T matrix. The numbers 
quoted for the pole position correspond to the pole in the u channel. The errors introduced 
by the expansion are, however, not relevant for the following discussion. In the commonly 
used notation {37], where the four sheets of the energy Riemann surface are defined as 

• Impi > 0, lmp 2 > 0: sheet I (physical); 

• Impi < 0, Imp 2 > 0: sheet II; 

• Impi < 0, lmp 2 < 0: sheet III; 

• Impi > 0, Imp2 < 0: sheet IV, 

this pole is located on sheet IV. A pole of this kind, being an analogue of a single chan- 

3611 . It should 



nel virtual state, was termed a "shadow resonance" by Eden and Taylor 
be distinguished from a "regular" Breit-Wigner resonance located on sheet III. The most 
important difference between a shadow resonance and a Breit-Wigner resonance is that, 
while the latter is accessible immediately from the physical region, i.e. the upper edge of 
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FIG. 2: Poles and zeros of the S matrix in the energy plane. Dotted line — a path from a point 
near the physical region, showed by a dot, to sheet III; dashed line — a path from the physical 
region to sheet IV. Roman numbers in brackets along the paths indicate the corresponding change 
of the sheets. The location of the pole in sheet IV, coinciding with the location of a zero of S22 on 
the physical sheet, is denoted by the circled cross. The location of the conjugated zero of £22 on 
the physical sheet is shown by the circled dot. 

the cut on the physical sheet, by moving down across the cut at KeE > A, the former can 
be accessed from the physical region only by moving around the zero-energy threshold, and 
is thus rather far away from the physical region. These different paths, corresponding to 
poles on sheets III and IV, are illustrated in Fig. |2j A very well-known consequence of a 
Breit-Wigner pole, lying on sheet III close to the physical region, is the expansion of the S 
matrix close to the pole, having the form 

iA 



S = I 

E-Er + ir/2' 

where A, the residue of the S matrix, is a Hermitian rank one matrix, satisfying 



(48) 



(49) 



where the quantities Tj are identified with the partial widths corresponding to the i th channel. 
These properties of A are dictated by the unitarity of the S matrix in the physical region. 
A pole on sheet IV, on the contrary, does not allow for a similar expansion of the S matrix 
in the physical region, due to the large distance to the pole. As a result, the residue, A, is 
no longer constrained to be Hermitian. The condition of Eq. (149]) is broken in this case as 
well: although the quantities \Au\ = r, can still be identified with the partial widths (see, 
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381]). they do not add up to the total width. In our case T p = 2.47 MeV, 



e.g., Refs. 128LI29 

T n = 0.34 MeV, corresponding to the so-called strength of the resonance ^jTj/T = 22.89. 

It is interesting to note that the pole, despite being located far away from the physical 
region, still crucially affects the S matrix and the observables close to the neutron threshold, 



via the zeros of certain elements of the S matrix and unitarity (see Ref. |c 



for a discussion 



of an analogous situation in the coupled system d(t,n) 4 B.e). Namely, associated with the 



pole on sheet IV, there is a zero of one of the e 
same energy on sheet I, as discussed in Ref. 



ements of the S matrix, S22, located at the 
j. Although this zero is on sheet I, it is still 



located rather far away from the physical region. However, it can be easily shown that there 
is another zero of S22, located at the complex-conjugate point, E = E r + ir/2 on sheet I. 
This zero is, in contrast to the pole and the zero discussed above, located very close to the 
physical region, as is also illustrated in Fig. |2j The nearby zero forces S22 (and, through 
unitarity, also S\\ above the neutron threshold) be small at E ~ E r , while the off-diagonal 
element, S12, approaches the unitary limit. These features are seen in the observables. 



IV. DISCUSSION 



In this paper, we have used an RG analysis to elucidate the power counting needed to 
organise an EFT describing a system of coupled channels with a single low-energy state 
close to one the thresholds. Such a system has, as its low-energy scales, the momenta 
corresponding to the various thresholds, as well as the on-shell momentum. These provide 
the expansion parameters of the corresponding EFT. 

The RG analysis of the two-channel case leads to the identification of three fixed points. 
One of these is just the trivial one. The expansion around it can be used to analyse systems 
that are weakly interacting in all channels. A second has two low-energy bound or virtual 
states and so both channels are strongly interacting at low energies. Lastly, of most interest 
in practice is one with a single bound or virtual state. Here, one linear combination of the 
two asymptotic channels is strongly interacting, the other weakly so. 

The power counting for the strongly interacting channel near that third fixed point is 
that of an effective-range expansion. As in the corresponding single-channel case, terms in 
the effective potential are promoted by two orders relative to simple dimensional analysis. 
The leading, energy-independent term generates large scattering length. It is a relevant 
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perturbation and so must be treated to all orders. All other terms in the potential can be 
treated as perturbations. 

All terms in the weakly interacting channel are of natural strength and so can be organ- 
ised according to simple dimensional analysis. The off-diagonal interactions that couple the 
two channels are promoted by one order. The leading one of these is then a marginal per- 
turbation, but it can be absorbed into the mixing angle that defines the strongly interacting 
channel. 

The resulting expansion can be applied to "halo" EFTs that describe nuclear systems 
with weakly bound low-energy states 5^^]. It may also be applicable to some of the states 
seen close to thresholds in the quarkonium systems 201 ] . 

In general the channels in these examples involve charged particles and so we have ex- 
tended our RG analysis to include Coulomb potentials between the particles. As in similar 
analyses of the single- channel case, these potentials are not sufficiently singular to substan- 
tially change the power counting. One obvious effect is to provide additional low-energy 
scales: the inverse Bohr radii for the various channels. 

In the strongly interacting channel, Coulomb forces introduce a marginal interaction. 
As a result, there is no longer a true fixed point but instead a logarithmically evolving 
RG trajectory. Theories near this trajectory can still be organised according to a modi- 
fied effective-range expansion in their strongly interacting channel. The weakly interacting 
channel is still natural, while off-diagonal interactions are again promoted by one order. 

To demonstrate the use of the resulting EFT, we have applied it to an ideal test case: 
the coupled p + 7 Li and n + 7 Be channels. They couple to a J p = 2~ state of 8 Be which lies 
within a few keV of the n + 7 Be threshold. 

At NLO we get a good description of the available p + 7 Li phase shift and the 7 Be(n,p) 7 Li 
reaction cross section using four parameters. At NNLO there is one further parameter, and 
although we are able to get fits of similar quality, we find evidence that the available data are 
not sufficient to determine all of the parameters at this order. The sizes of the perturbative 
terms are consistent with an underlying scale of the order of 50 — 100 MeV. The scattering 
length in the strongly interacting channel is about 0.1 MeV -1 , making it unnaturally large 
on this scale. The differences between the paramaters from fits at LO, NLO and NNLO are 
also as expected. 

As in previous analyses of this system, we find that the J p = 2~ state of 8 Be is described 
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by a pole on the fourth sheet. This lies far from the physical sheet so it is not a resonance but 
rather the multi-channel analogue of a virtual state. It makes its presence known through 
the cusp in the p + 7 Li phase-shift at the n + 7 Be threshold and the very large cross section 
for 7 Be(n,p) 7 Li at low energies. 

The example demonstrates the viability of this EFT approach for analysing coupled- 
channel systems with low-energy bound or virtual states. It may provide a more systematic 
alternative to the i?-matrix method which is widely used in studies of such systems. It can 
also be applied to similar states in the quarkonium systems. 
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